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Phase contrast imaging is used to observe Bose-Einstein condensates (BECs) at finite temperature 
in situ. The imaging technique is used to accurately derive the absolute phase shift of a probe 
laser beam due to both the condensate and the thermal cloud. The accuracy of the method is 
enhanced by using the periodicity of the intensity signal as a function of the accumulated phase. 
The measured density profiles can be described using a two relevant parameter fit, in which only the 
chemical potential and the temperature are to be determined. This allows us to directly compare 
the measured density profiles to different mean-field models in which the interaction between the 
condensed and thermal atoms is taken into account to various degrees. 



PACS numbers: 03.75.Hh, 67.85. Be, 87.64.mh 

I. INTRODUCTION 

In the field of ultra-cold atomic physics the two com- 
monly used techniques to image a cloud of atoms are 
absorption imaging and fluorescence imaging, i. e. imag- 
ing the cloud's absorbed or radiated intensity, respec- 
tively. In absorption imaging the shadow cast in a probe 
beam by the cloud of atoms is recorded on a CCD cam- 
era. From the resulting image the spatial extent of the 
cloud and the optical density, proportional to the col- 
umn density along the line-of-sight of the probe beam, 
can be derived. For dense clouds absorption imaging be- 
comes unreliable, since for optical densities of the order 
of 4 the dynamic range of the absorption imagining tech- 
nique is insufficient to accurately determine the intensity 
in the shadowed region [T]. For Bose-Einstein conden- 
sates (BECs) the typical in situ optical density is of the 
order of 500 and the image gets completely "blacked out" , 
making it impossible to extract the column density from 
the image. The fluorescence imaging technique critically 
depends on the intensity and the frequency of the probe 
beam. Furthermore, for high densities the cloud becomes 
opaque and atoms in the center of the cloud are not sub- 
jected to the same intensity as atoms on the edge of the 
cloud. 

One way to reduce the optical density is to reduce 
the atom-photon cross section by choosing a large probe 
beam detuning. However, for nonzero detuning the real 
part of the index of refraction of the atoms becomes 
nonzero and the cloud behaves like a gradient-index lens. 
This results in a refraction of the probe beam and leads 
to distortion of the images. The resulting images cannot 
be used for a relative or absolute column density mea- 
surement, since the magnitude of this effect has a spatial 
dependence. Detuning the probe beam further from reso- 
nance reduces the index of refraction but leaves the cloud 
virtually transparent before the refraction is sufficiently 
reduced. 

The optical density can successfully be decreased by 
reducing the atomic density. The cloud is allowed to ex- 
pand during a specific time-of-flight (tof ) and an image is 



recorded as soon as the optical density is of the order of 2. 
Although absorption imaging after expansion is the com- 
monly used technique in this field, expansion of the cloud 
introduces some difficulties as well. First, if the cloud is 
a BEC at finite temperature there is no accurate descrip- 
tion available of the expansion of the cloud. Second, the 
quantization axis of the cloud is no longer well-defined 
since the magnetic fields are turned off. As a result, 
the effective cross section for the polarized probe beam 
becomes ambiguous. As a consequence both absorption 
imaging and fluorescence imaging cannot be applied to 
accurately determine the column density of dense clouds. 

Some techniques have been developed to image dense 
clouds by using strong saturation in both absorption [2] 
and fluorescence [3] imaging. However, the high-intensity 
regime has to be calibrated based on the low-intensity 
regime and the latter regime is imaged after expansion. 
As a consequence the issues raised above about imaging 
after expansion apply to the strong saturation techniques 
as well, although cloud can be imaged in situ. 

Alternatively, phase contrast imaging (PCI) has been 
successfully applied to image dense atomic clouds [H |5J 
6 . In PCI the real part of the index of refraction is 
measured. The detuning can be chosen in such a way 
that the refraction is small but the phase shift due to 
the index of refraction remains substantial. PCI can be 
used to image BECs in a nondestructive way, since the 
number of scattered photons becomes negligible for large 
detuning. Other techniques based on dispersive light 
scattering, such as diffraction-contrast imaging of cold 
atoms [7] and nondestructive spatial heterodyne imag- 
ing of cold atoms [S] have been demonstrated on dense 
thermal clouds. 

Here, we use PCI to measure the in situ density dis- 
tribution of dense BECs at finite temperature absolutely 
and accurately. By choosing a smaller detuning com- 
pared to the nondestructive PCI schemes, the density 
of the thermal cloud and the BEC are determined ac- 
curately for all typical temperatures. This allows us to 
directly study the interplay between the condensed and 
thermal atoms and the influence of this interaction on 
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the transition temperature for condensation. 

This paper begins with a detailed description of the 
PCI technique and the way it is appHed to a cold Bose 
gas. Next, the accuracy is demonstrated in two experi- 



ments. In the first experiment, described in § VCl the 
number of condensed atoms is measured with an accu- 
racy below 5% in a single shot. In the second experiment, 
described in § |VD| the growth of the condensate is stud- 
ied as a function of the temperature. In this measure- 
ment , the interaction between the condensed and thermal 
atoms is studied directly. Furthermore, measurements 
around the temperature where a BEC is formed are used 
to determine the transition temperature. 



II. THEORETICAL DESCRIPTION OF PCI 
A. Imaging the phase 

It is possible to image a cloud of atoms by illuminat- 
ing the cloud with a probe beam and record the shadow 
cast due to the absorption of the atoms on a CCD cam- 
era. Since the optical density is measured, this technique 
probes the imaginary part of the complex index of re- 
fraction J\f of the atomic cloud. In contrast, the imaging 
technique described in this paper probes the real part of 
Af. This technique is used in the phase contrast micro- 
scope which is developed by Zernike in 1933. In 1953 
Zernike was awarded the Nobel prize "for his demon- 
stration of the phase contrast method, especially for his 
invention of the phase contrast microscope" ^ . 

In general, the complex phase </> accumulated by the 
light in a medium such as a cloud of atoms can be ex- 
pressed as 0atoms = 4> + i4> /2, with 4> the optical 
density. Writing the electric field of the probe beam as 
Eprohc — Eq exp{i(f)) with Eq the amplitude of the light 
field, the electric field after passing a cloud of atoms can 
be written as 
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where the electric field is split in two parts: the part that 
is diffracted by the atoms, -Eatoms, and the part that is 
not, iJprobc- In absorption imaging the phase information 
is lost on a CCD camera, since the intensity / is mea- 
sured as / ceoEE*/2 = /oe"'^" where !„ = ceo\Ea\'^/2, 
c is the speed of light and Eq is the vacuum permittiv- 
ity. So the real phase term (/) cancels and only the ab- 
sorption is measured. In the simplest form of PCI, the 
non-diffracted part of the probe light (first term on the 
RHS of Eq. 0) is blocked. The light that arrives at 
the camera is the light which is diffracted by the atoms, 
E = -Eprobo(exp (z(/)atoms) " !)• Now the phase informa- 
tion no longer cancels and the intensity on the camera is 
given by 



/ = 2/o(l - cos (0atoms)), 



(2) 



where the arbitrary phase of the probe beam cancels. 
This technique is referred to as dark-field imaging and 



has been successfully used to image ^^Na atoms in the 
early days of experimental BEC physics [6] . An impor- 
tant experimental drawback of this method is the absence 
of probe light on the camera preventing the normaliza- 
tion of the intensity profile, which makes this technique 
not suitable to measure the absolute phase shift. How- 
ever, the relative signal can be used to study the spatial 
distribution. 

A more elaborate way of converting the accumulated 
phase into an intensity profile is by phase shifting the 
non-diffracted light instead of blocking it, in the same 
way as the phase contrast microscope works. The non- 
diffracted light is the plane wave part of the beam, which 
can be phase-shifted by placing a small transparent ob- 
ject with a refractive index in the probe beam in such a 
way that only the plane wave part of the beam propa- 
gates through it. This object will be called a phase spot. 
The plane wave part propagating through the phase spot 
accumulates a phase 6, changing Eq. ([ij to 



Now, the intensity is given by 



!)■ 



(3) 



/= /o(3-2cOs(6')-|-2cOs(e'-(/>atoms)-2cOs(0atoms)), (4) 

where the arbitrary phase of the probe beam cancels. The 
probe beam is not blocked and it is possible to derive the 
absolute phase shift, since the intensity profile Iq can be 
measured and used for normalization. This technique has 
also been successfully applied to image ^'^Na atoms [S]. 
The implementation of PCI described in Ref. [5] uses a 
detuning of the probe beam, which is large compared to 
the atomic resonance and the absorption of the light by 
the atoms is therefore reduced to a degree that it is no 
longer destructive. This allows for taking multiple images 
of the same BEC. 

In our experiment we do not focus on the ability to 
make a nondestructive image but we use the periodicity 
of the intensity as a function of (/fatoms (see Eq. Q) to 
accurately determine the accumulated phase of both the 
BEC and the thermal cloud. In recent years, different 
schemes are proposed derive the phase information of a 
cloud of cold atoms [SJ |TU] , but the accumulated phase 
has never been used for accurate quantitative measure- 
ments of the density distribution of BECs at finite tem- 
perature. 



B. Index of refraction of a cold gas 

We continue with a detailed description of the phase 
shift of an electro-magnetic wave propagating through a 
cloud of spin-polarized atoms in order to relate (/)atoms 
to the atomic density. The phase shift is caused by the 
modification of the probe field by the atomic dipoles, so 
we start with an expression of the polarizability of one 
atom. The polarizability tensor in the low-intensity 
limit and (5 ^ w is given by |llj 
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(g|/Zgo|e)(e|^cg|5')(5Vggl5), (5) 



where the natural hnewidth 7 is given by 



(6) 



Here, fieg,ge are the electric-dipole operators, Cgg is the 
density matrix, 5 is the detuning and u is the probe fre- 
quency. In the experiment various Hnear polarizations of 
the light field have to be dealt with, so the dependence of 
the polarizability ^ on the polarization has to be deter- 
mined. We define the angle (3 as the angle between the 
z-axis of the atomic cloud and the (linear) polarization 
direction of the light field which is in the (x, z) plane. In 
this system, the line-of-sight is along the y-axis. In order 
to evaluate Eq. ([5| for different polarizations, a fixed po- 
larization axis is chosen parallel to the polarization of the 
light field and Eq. ([5| is evaluated in this frame, in which 
conveniently only AM — transitions are induced. Fur- 
thermore, only one component of the polarizability ten- 
sor will contribute in this frame and the polarizability 
becomes a scalar. Since all atoms in the trap are in the 
|Fg, Afg) = |1, —1) state with respect to the z-axis of the 
magnetic field confining the atoms, the density matrix 
Ugg only has one nonzero element, which makes the cal- 
culation of the rotation from the magnetic z-axis to the 
polarization axis of the light field straight-forward. 

Rotating the density matrix of the ground state agg is 
achieved using 



,(/3) = 7^,(/3)tfT3,7^,(/?), 



(7) 



where the rotation matrix TZy(P) for rotating a system 
around the y-axis over an angle (3 for a system with Jg — 
1 can be found using Wigner's formula as [12j 
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Since the two dipole moment operators fleg do not 
change the M-state of the atoms, only the diagonal ele- 
ments of Eq. ([5]) play a role and we end up with 
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2t ^ (Fg,Mg|^ge|fc,Mc)^(-Fg,A/g|»gg(/3)|Fg,Mg 
h^j^ 1 - 2z5e/7 



(9) 

Furthermore, we can define the square of the Clebsch- 
Gordan coefficients Cp^.M^ (transition strength) as 



and obtain using Eq. ^ 

^ _ leocax ^ CgAF^.M^\aggmF^,M^) 
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where 5e is the detuning with respect to transition Fg 
Fc and a\ the cross section for absorption of light: 
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For sodium atoms in the |i^g,Mg) = |1,~1) state this 
results in 



where 
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(13) 

(14a) 
(14b) 
(14c) 



For alkali metal atoms the polarizability is independent 
of the angle /3 if the detuning S is large compared to the 
hyperfine splitting of the excited state. In that limit the 
numerator of the sum in Eq. ( 13 1 is | with no angular 
dependence. Finally, the complex mdex of refraction Af is 
given by 



AA2 = 1 



pa 
£0 



(15) 



with p the density. Note that for atoms confined in a 
potential the density is not homogeneous and p depends 
on the position, p = p{x, y, z). 



The index of refraction given by Eq. ( 15 ) is only valid 



for low densities. For higher densities, the dipole moment 
of an atoms is influenced by the internal field of the sur- 
rounding atoms. This results in a modification of the 
refractive index, which is accounted for by the Lorentz- 
Lorenz equation |13j 



where C is given by 



C = --pa/eo. 



(16) 



(17) 



Since the atoms being imaged are bosons instead of 
classical particles, the index of refraction will be modi- 
fied further due to bunching of the atoms. Two terms 
are added to Eq. (17). The first term accounts for mod- 



ifications due to the enhanced photon scattering cross 
section caused by the quantum statistics of the atoms. 
The second term describes the modified refractive index 
induced by resonant Van der Waals interaction, resulting 
in an increase of almost 10% in both the real part and 
the imaginary part of the refractive index |14j . 

We have derived the modifications along the lines of 
Ref. [14j and have confirmed the result for the detuning 
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and densities reported in that paper. We found modifica- 
tions to the complex index of refraction strongly depend 
on the chosen detuning and therefore on a. In Ref. [H] . 
i5 ~ 7, while we are in the regime 5^7. For our param- 
eters, the real part of the index of refraction does not 
change significantly. However, we do find an enhance- 
ment of the absorption by a factor of 3. Since the overall 
absorption is very small, up to four percent at the detun- 
ing used, enhanced absorption does not influences the 
intensity profile significantly. As a result, we expect the 
measured phase to change less than one percent due to 
this effect. 



C. Imaging the atomic density 

The image of a cloud of atoms recorded on a CCD 
camera yields the column density of the cloud, that is the 
density profile integrated along the probe linc-of-sight. 
The column density can be linked to the complex index of 
refraction in the following way. In general, light passing 
through a medium with an index of refraction J\f over a 
distance I accumulates a complex phase ((> given by 



(18) 



where k — 2t:/X. Since the index of refraction depends 
on the position, the accumulated phase is written as 




p{x,y, z)a/eQ 
l + C 



1 dy. (19) 



Since jpa/eol *C 1 the index of refraction can be approx- 
imated as 



pa/ So ^ pa 
l + C ~ 2eo ■ 



(20) 



Under typical experimental conditions |pa/eo| is 10~^, 
introducing only a small error by using this approxima- 
tion. Now, Eq. ( 19 1 can be written as 



a 1 a 

(t){x,z) = k— / p{x,y,z)dy ^ k—pc{x,z), (21) 
2eo J 2eo 

where the integration is along the linc-of-sight y and 
Pc{x, z) is the column density. 

The real part of (/)(a;, z), proportional to Re (a), yields 
the phase shift, whereas the imaginary part, proportional 
to Im(a), yields the absorption. Since a is constant 
for fixed detuning 5, the phase shift is directly propor- 
tional to the column density. Since Eq. ( 13 ) satisfies the 
Kramers-Kronig relation 



Re (A/") - 1 



Im (A/") , 



(22) 



the ratio between the phase shift and the absorption is 
constant for a given detuning. 
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FIG. 1; Schematic representation of the PCI setup. The 
atoms located at the center of the trap are imaged with lens 
Li with focal distance /i which is placed at a distance /i from 
the center. A sharp image is created by placing the lens L2 
at distance d from L\ with focal length /2 at a distance /2 
from the CCD camera. The phase spot (PS) is placed in the 
focal plain of the non-diffracted probe beam. 



III. PCI SETUP 

In order to create a phase contrast image of the 
atoms, the plane wave part of the probe field has to be 
phase shifted as is described in § |II A[ This is commonly 
achieved by Fourier transformation using a lens. The 
concept of the imaging setup used in our experiment is 
schematically shown in Fig. [T] A two-lens setup is used 
to generate a sharp image of the atoms on a camera. The 
first lens Li with focal distance /i — 250mm is placed at a 
distance /i from the center of the magnetic trap. The sec- 
ond lens L2 with focal distance J2 = 750mm is placed at a 
distance /2 from the CCD camera. The distance between 
both lenses is d = 150mm. The combination of both 
lenses creates a sharp image of the cloud on the CCD 
camera with a magnification M — /2//1 = 3.0 ± 0.05 
where the uncertainty is the result of the uncertainty in 
the focal distance of the lenses. 

In order to shift the phase of the plane wave part of 
the light a phase spot is placed in the back focal plane 
(Fourier plane) of the two lenses, which lies for these 
parameters beyond the second lens L2. This position 
can be found easily in the experiment, since the parallel 
probe beam has its waist at the back focal plane and 
therefore the phase spot is placed at this position. 

Figure [2] shows the intensity I / Iq as a function of the 
accumulated phase ^atoms for different values of 9 based 
on Eq. Q. In all experiments described here, = 7r/3 
is used. This value is chosen since it yields a maximum 
visibility (/max - /min)/ (/max + /mill) = 1 and a dynamic 
range /max ~ ^min = 4/9. A visibility of 1 is reached only 
for 6 = 7r/3, as can be seen in Fig. [2) Further more the 
response is approximately linear for small accumulated 
phases. 

For technical reasons, the 7r/3 spot consists of a glass 
plate with a dimple instead of a spot on top. This dim- 
ple is made by dry-etching a 50fim diameter round dim- 
ple with a depth of 5A/(6(7V - 1)) = 1071 ± 15nm in 
a 50mm round fused silica glass plate of 4mm thickness 
\15\ . Light propagating through the dimple accumulates 
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FIG. 2: The phase contrast signal I/Io as a function of 
the accumulated phase (/>atoms for different phases 9 = 
TT, 7r/2, 7r/3, 7r/4 of the phase spot. 



57r/3 less phase than the light propagating through the 
full thickness of the plate and thus 6 = +7r/3. 

The atoms are imaged using an Apogee APIE camera 
with a Kodak KAF-0401E chip camera with a pixel size 
of 9 X 9 jjLm? . At the magnification M = 3.0, the effective 
camera resolution of 3.0 /im per pixel is comparable to 
the diffraction limit Xros ~ l-22A/i/(2r) « 3.6 /xm of the 
imaging lens Li with radius r = 25mm for probe light 
with a wavelength A = 589nm. 



IV. DENSITY DISTRIBUTION 
A. Three models for BECs at finite temperature 

A phase contrast image of a BEC at nonzero tempera- 
ture contains the column density distribution of the two- 
component cloud, i. e. a BEC and a thermal cloud which 
can be distinguished from each other due to their distinct 
density distributions. In order to derive the relevant pa- 
rameters from such an image, like the thermal (excited) 
density ricx, the condensate (ground state) density ric, 
and the temperature T, the measured profiles are fitted 
to a theoretical model describing the cloud. Various mod- 
els exist in literature in which the interaction between the 
condensed and thermal atoms are taken into account to 
various degrees. In this paper we compare our measure- 
ments to three of these models. The first model is the 
commonly used bimodal distribution in which all inter- 
actions between the two components are ignored. This 
model is used in to describe the vast majority of the ex- 
periments of cold Bose gases In the second model 
only the effect of the condensed atoms on the thermal 
atoms is taken into account. The third model incorpo- 
rates the effect of the density distribution of the thermal 
atoms on the condensed atoms as well. Below, the den- 
sity distribution of the thermal cloud and the BEC in 
each of these models are described consecutively. 

In the first model, referred to as the ideal model, the 
chemical potential is given by /i = Vc^t{r) + nc{r)Uo, 



where Uq = Anh^a/m is the interaction parameter for 
bosons with s-wave scattering length a and mass m |16j . 
Since there are no interactions between condensed atoms 
with density Uc and thermal atoms with density rigx, M 
does not depend on riex- In the ideal model riox is given 

by 



where e is given by 



dp 



1 



(27r;i)3e('^W-^)/'=B7^- 1' 



2m 



(23) 



(24) 



In this approximation both riox and ric can be derived an- 
alytically yielding the bimodal density distribution, the 
sum of a Maxwell-Bose distribution modeling the thermal 
cloud and a Thomas-Fermi (TF) distribution modeling 
the BEC. In this approximation nox is given by 



Li3/2Ue-^-'(^^)/'=-^)/A: 
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(25) 



where the de Broglie wavelength is given by AdB = 
■y/ {2TTh?) I (mk^T) and z — exp{fj, / (kBT)) is the fugac- 
ity. Li3/2(2) is the Bose function, given by Li„(z) = 
Er=i ^Vfc" and 



/ 2 2 I 2 2 I 2 2\ 



(26) 



is the external potential. In the TF approximation ric is 
given by 



Uo 



■Tx 
Riad 
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(27) 

where i?rad ax are the TF radii of the condensate given 

by 



-Rrad — 



-Rax — 



2^ 



rad 
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(28a) 
(28b) 



Since Wrad/^ax ~ 50 in the experiments described here 
and thus i?rad ^ -Rax the clouds are cigar shaped. The 
number of condensed atoms Nc is found by integrating 
Tie over f and yields 



N, = / n,{r)dr = 



Sir f 2fi 



15 



3/2 



(29) 



Integration of the density over the line-of-sight yields the 
column density, which is the property being measured in 
the experiments. 

The second model, denoted as the semi-ideal model, 
does include interactions, but only the contribution of 
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the mean-field potential of the condensate to the thermal 
density distribution is taken into account [17 . Therefore, 
fi remains unchanged: /i = Vc^ti^ + "■c(^C^i but the 
energy of the thermal atoms is changed due to the mean- 
field potential of the condensed atoms and written as 

e^{r) = (pV(2m) -I- 2n,{r)Uo + T/ext(r) - fi)^ ~{n,{r)Uof 

(30) 

We calculate the density distribution in this semi-ideal 
model by numerical integration of Eq. (23), where e{r) 
is given by Eq. (30) [16]. The shape of the density dis- 



tribution of the condensate remains unchanged, but the 
repulsion of the mean-field potential of the condensate 
causes the thermal atoms to be repelled from the center, 
which results in the shape shown in Fig. [3] In order to 
have the same total number of atoms and temperature 
in Fig. [3j the chemical potential is different in both mod- 
els. The decrease in density in the center is the result of 
the repulsion due to the mean-field potential of the con- 
densed atoms and therefore absent in the noninteracting 
distribution. Due to the lower central density the result- 
ing condensate fraction is lower in the semi-ideal model 
compared to the noninteracting model for the same total 
number of atoms N — A'cx + -^c ^^nd temperature T. 

In the third model the mean-field effect of the ther- 
mal cloud on the condensate is taken into account as 
well. This model corresponds to the Popov approxima- 
tion in the TF limit, so we refer to this model as the 
Popov model The chemical potential is now given 
by /z = Vext{r) + [ndf) + 2nox(?^] Uq- The expression for 
the density of the excited atoms in the Popov approxi- 
mation is given by 



dp pV2m + 2n(r)[/o + Kxt(r)-M 



Qe(r)/kBT 

(31) 

with n = Hex + nc and where, in the semi-classical ap- 
proximation, e(r) is given by 

e\r) = (//(2m) + 2n{r)Uo + V,xt{r^ - fi)^ -{n,{r)U„y 



(32) 

Now, the resulting density distribution of both ricx and 



Uc is found by numerical integration of Eq. (31 1, where 



an iterative procedure is used to find a self-consistent 
solution. The difference between the semi-ideal model 
and the Popov model is shown in Fig. |4] focusing on the 
density distribution of the condensed atoms. The den- 
sity distribution of the thermal atoms as shown in Fig. [3] 
results in a higher effective potential for the condensed 
atoms. This effect results in a compression of the density 
distribution of the condensate compared to the distribu- 
tion calculated using the semi-ideal model. 



B. Determining the chemical potential in the 
experiment 

In the theoretical models the density profile is deter- 
mined by the parameters ^ and T in a harmonic con- 




Position (mm) 



FIG. 3: The density distribution of a thermal cloud Wex cal- 
culated using the noninteracting model (dashed line) and the 
semi- ideal model (solid line), for T = 0.9/iK, A'^ = 1- lO^atoms 
for the typical experimental trap parameters. The inset shows 
the total density distribution ricx + nc for both models. 
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FIG. 4: The combined density distribution of the thermal 
cloud and condensate riex -|- fic calculated using the semi-ideal 
model (dashed line) and the Popov model (solid line) for T = 
0.9/ jK. A'' = 1 • lO^atoms and the typical experimental trap 
-paiameters. The left and right insets show an enlargement of 
the edges and the center of the condensate, respectively. 



finement given by Vcxt- Therefore, it seems evident to 
analyze the measured density profiles by fitting them to 
a theoretical model in which only /i and T are to be de- 
termined. However, this approach can be applied only 
if the size and the height of the measured profiles can 
be determined absolutely. Since the aforementioned PCI 
method is expected to do so, we can apply such a fit 
procedure here. 

In the commonly used imaging methods the size or the 
height of the density profiles is not measured absolutely, 
and one has to fit the size and height of signal of both 
the BEC and the thermal cloud separately. In that case, 
fi and T are derived from the height and size of the mea- 
sured density profiles. The experimental conditions favor 
the use of one of these measures above the others. 

Being able to describe the complete cloud by fi and T, 
as is done in the theoretical models, is expected to yield 
a higher accuracy, since different measures for the chem- 
ical potential and temperature are combined which de- 
pend differently on experimental parameters such as the 
magnification, detuning and the trap frequencies. Fur- 
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thermore, it allows for taking the interaction between 
the condensed and thermal atoms into account. This in- 
teraction is ignored in the vast majority of experiments 
in this field, in contrast to the theoretical description of 
a condensate at finite temperature found in many text- 
books. Our PCI setup is expected to yield an accurate 
and absolute determination of both the density and the 
size of the cloud and measuring in situ prevents the need 
to account for the expansion of the cloud. Therefore, 
since all measures are absolute, fi and T are expected 
to fully describe the measured density profiles, given the 
fixed experimental parameters such as the magnification 
and the trap frequencies. This allows us to directly com- 
pare our measured density profiles to the three mean-field 
models described in the previous section. 

The PCI method not only yields absolute density pro- 
files, it is expected to yield a more accurate determination 
of the thermal fraction for low temperatures compared to 
absorption imaging as well. Since the accumulated phase 
for the dense BEC can be multiple times tt, a significant 
phase is accumulated in the thermal cloud, without sat- 
urating the signal in the BEC. In absorption imaging on 
the other hand, the optical density of the thermal cloud 
is strongly reduced in order to reduce the optical density 
of the BEC to the order of 4. Especially for low tem- 
peratures, when only the low density tail of the thermal 
distribution spatially extends the BEC, it turns out to be 
difficult to distinguish between both components in the 
experiment. Various schemes have been used to enhance 
the contrast between the thermal cloud and the conden- 
sate, for example by making a two-pass fit routine, first 
only on the thermal cloud, followed by the condensate 
and thermal cloud while fixing the temperature found in 
the first pass [TS]. Another scheme used is spatial sepa- 
ration of both components using Bragg spectroscopy and 
fitting both parts separately [TB]. In the PCI method 
the contrast between both components is sufficient for 
all measured temperatures and prevents the need to in- 
troduce elaborate schemes to enhance the contrast. 



V. EXPERIMENTAL RESULTS 

A. Experimental parameters 

Imaging the atoms is conducted in such a way that 
the periodicity of the intensity of the PCI technique is 
used. Since the intensity signal varies periodically as a 
function of (/>atoms (Eq- Q), the intensity signal shows 
rings in the intensity profile for sufficiently large values 
of 0atoms (see Fig. |2]). The number of rings depends on 
(^atoms, which scales with the density and the detuning 
S as can be seen from Eqs. (211 and 13 



The detuning 

is a parameter which can be tuned accurately in a wide 
range (up to 400MHz) using acoustic-optical modulators 
(AOMs) allowing us to make images with an adjustable 
number of rings. Due to the ring pattern the imaging 
lens can be put in position very sensitively allowing us 



to resolve features down to 4/im. Furthermore, S can 
be chosen in such a way that the intensity profile of the 
BEC shows rings and the less dense thermal cloud yields 
a significant intensity signal as well, even for low tem- 
peratures. The resulting effective dynamic range of this 
method is therefore increased by the periodic dependence 
of the intensity on the accumulated phase. 

A drawback of this method lies in the fact that the 
BEC acts as a gradient index lens, the lensing being 
stronger for smaller 6. The effect of the lensing on the 
imaging resolution is estimated by a computer simulation 
in which the Fresnel-Kirchhoff diffraction integral is nu- 
merically integrated. In this simulation the finite size of 
the imaging lenses is explicitly taken into account. The 
mutual distances between the BEC, lenses, phase spot 
and image plane are chosen identical to the values used 
in the experiment. The intensity distribution at the im- 
age plane is used to estimate the effects of the refraction 
of the BEC on the imaging resolution. If the diffraction 
is too large, the higher diffraction orders can miss the 
aperture of the imaging lenses and will degrade the reso- 
lution of the image. This simulation yields an upper limit 
for the allowed phase shift of the BEC and therefore a 
lower limit of the applied detuning given the condensate 
density and dimensions. In the experiment we are in the 
regime in which the aperture of the lenses is not limiting 
the resolution, although the onset of a slight deformation 
of the density profile is observed for the highest conden- 
sate densities. For the typical parameters, the maximum 
focal length of the BEC is found to be of the order of 
100/iTO. Even though the diffraction is not too large for 
the apertures in the imaging path, the imaging lenses 
have to be placed in focus within a few tens of microme- 
ters to prevent distortions of the imaged intensity profile, 
since the imaged object itself acts like a lens. 

The phase spot has to be aligned accurately as well, 
since all of the plane wave light has to propagate through 
the phase spot. In the long axis of the cloud the diffrac- 
tion is weak, and for too large a phase spot the diffracted 
light propagates through the phase spot as well. In order 
to make a off between the two criteria that the phase spot 
should be large enough that the non-diffracted light is 
phase-shifted, but small enough that the diffracted light 
is not, the phase spot has a size of only a few times the 
waist of the probe beam. As a consequence, the phase 
spot is easily misaligned. Misalignment of the phase spot 
leads to reflection of the probe beam on the edges of 
the phase spot and possibly shifting the phase of the 
light, which is diffracted by the atoms. As a result of the 
misalignment we notice a higher, tilted intensity profile, 
which is easily detected, since we have several maxima in 
the intensity profile. The common way to apply PCI is 
by using a large detuning and the intensity signal never 
reaches the first maximum value. Without a well-defined 
maximum in the image, (slightly) more intensity due to 
misalignments may be unnoticed and the method can no 
longer be used to determine the absolute phase. 

As can be seen in Fig. [2j the intensity as a function 
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of the phase varies in certain ranges rapidly. Choosing 
a detuning in such a way that the maximum phase lies 
in this range, a shght increase in the density leading to 
a slight increase in the maximum phase leads to a large 
change in the intensity in the center of the cloud making 
the method very sensitive to small density changes. 

In the experiment imaging the atoms consists of taking 
a sequence of three images. The first image obtains the 
intensity profile /atoms (a;, of the probe field and the 
atoms. The second image contains the intensity profile 
Io{x, z) of the probe field in absence of the atoms. This 
image is generally shot after the magnetic trap has been 
turned off for two seconds and no atoms are left. The 
third image is the intensity profile /bg(a;, z) with neither 
the probe field nor the atoms and is used as a background 
image. The normalized intensity profile /(a;, z) of the 
atoms is given by 

1{X,Z} — — r - - r , {66} 

Io[X,Z) - Ihg{X,Z) 

where the subtraction and division is performed on a 
pixel- by-pixel basis. The magnification scheme used 
causes the probe beam to get blown up to such an extend 
that the beam profile is no longer smooth on the cam- 
era due to interference effects caused by dust and small 
scratches on the imaging optics. In principle, these im- 
perfections are canceled in the final images, but since the 
imaging path is not interferometrically stable in between 
the acquisition of the three images, some distortions re- 
main. 



B. Phase contrast images 

To demonstrate the PCI method, atoms are cooled be- 
low the critical temperature in a 95.57Hz x 95.57IIz x 
2.234Hz trap and imaged with a probe beam with an in- 
tensity /o = 6-10^^ /sat for 50 lis. The probe beam is 
detuned 28. 17 below the |Fg = 1) ^ |Fc = 1) transition 
and results in typical intensity profile /(a;, z) as shown in 
Fig.|5] The high frequency noise found in these images is 
caused by the distorted profile of the probe beam. Small 
fringes close to the BEC are the result of the lenses being 
slightly out of focus. 



C. Accurate determination of the number of 
Bose-condensed atoms 

1. Accuracy using absorption imaging 

The usual way to determine the number of condensed 
atoms is by taking a series of absorption images in time- 
of-flight. Absorption imaging turns out to become un- 
reliable when the optical density exceeds 4 due to the 
limited dynamic range of the CCD camera [T]. If the 
cloud of atoms is a BEC, the typical optical density is in 



the order of 500 when the probe is on-resonance. Imag- 
ing the atoms off-resonance to take advantage of the re- 
duced photon scattering cross section turns out to be 
complicated since the cloud behaves as a gradient-index 
lens in this regime. The optical density can be lowered 
by turning off the confinement causing the atoms to ex- 
pand during a certain time-of-flight, until the density is 
low enough that the optical density is in the order of 4. 
However, expansion complicates the interpretation of the 
measured density profiles, since it consists of the convo- 
luted momentum and spatial distribution of the atoms 
and the expansion of the cloud cannot be described ex- 
actly. 

Furthermore, the switching of the magnetic confine- 
ment complicates the retrieval of the column density due 
to the finite time needed to switch the magnetic fields. 
In the absence of the magnetic field, the quantization 
axis of the atoms is no longer well defined causing the 
atoms to align along small residual magnetic fields dur- 
ing time-of-flight. This changes the effective cross sec- 
tion of the atoms for the applied polarized probe light. 
The magnitude and direction of the residual magnetic 
fields are expected to be spatially dependent. Since the 
atoms fall during time-of-flight due to gravity, the ef- 
fective cross section is also expected to depend on the 
time-of-flight duration. Furthermore, the imaging lenses 
have to be repositioned for each time-of-flight duration. 
A final complication is that the accuracy of detuning of 
the probe light is limited to IMHz, corresponding to a 
reduction of the absorption up to 4%. All these effects 
cause uncertainties in the measured number of atoms up 
to 20%. 

A technique is proposed based on the nonlinear re- 
sponse of atoms on the applied probe intensity in order 
to lower the optical density of the cloud without the use 
of expansion. This allows for quantitative in situ absorp- 
tion imaging, although the precise calibration required 
has to be done after time-of-flight, reducing the accuracy 
of this method [5]. Moreover, the required reduction of 
the optical density of the condensate results in a reduc- 
tion of the optical density of the thermal cloud as well, 
complicating the determination of both density profiles 
simultaneously. The high intensity needed to reach the 
strong saturation regime, / (x OD x /gat ~ 8W/cm^, with 
/sat the saturation intensity, complicates the implemen- 
tation of the method as well. 



2. Accuracy using PCI 

Using the PCI technique circumvents the issues raised 
above. The detuning is large compared to the uncer- 
tainty, 6 — 281 ± IMHz, and since the measurement is 
conducted in situ the confinement is not switched off and 
the quantization axis of the atoms remains well defined. 
The Zeeman shift caused by the confinement is less than 
2MIIz and therefore negligible compared to the detuning. 

An image of a BEC taken with PCI technique yields 
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FIG. 5: A set of four phase contrast images during the final stage of the evaporation process: (a) [ijh = 2.6kHz, T = l.Ol^K, 
iV = 7- lO^atoms, ISI^jN < 1%; (b) /i/h = 4.4kHz, T = 0.95/iK, iV = 6- 10** atoms, N^/N ^ 15%; (c) fi/h = 5.2kHz, T = 0.58^K, 
iV = 4 • lO^atoms, NJN « 50%; (d) fi/h = 4.8kHz, T = 0.25/iK, N = 2- lO^atoms, NJN > 95%. 



the accumulated phase and the size of the condensate. 
Analyzing all known systematic and statistical errors 
show that the largest contribution to the error in the size 
originates from the magnification of the imaging system. 
We determined the magnification to be M = 3.05 ±0.05. 
Since the condensate radius R is proportional to M, the 
number of condensed atoms scales as oc (see 



Eq. ( 28b I and Eq. ( 29 1) . The trap frequencies are derived 
from a center-of-mass oscillation measurements yielding 
the trap frequency with a statistical uncertainty below 
10-^ 

The largest contribution to the error in the accumu- 
lated phase originates from the lensing effect of the con- 
densate and is estimated to be less than 5% in the mea- 
sured phase based on the computer simulation we con- 
ducted. This effect does not influence the axial size of 
the cloud since the density varies slowly in this direction 
making the lensing negligible. Since the number of atoms 



scales with the accumulated phase as cx cf)^^^, lensing 
is expected to yield an error up to 8% in the number of 
atoms for the highest densities. 

We analyze measured clouds for similar parameters as 
the ones shown in Fig. [5] where the condensate fraction 
is at least 90% by making a least square fit to the ideal 
model. We conclude both the determined size and phase 
of the condensate yield the same chemical potential fi 
within 3% for the typical number of condensed atoms 
Ac ^ 2.5 • 10^. Since the thermal fraction is small, inter- 
actions are expected to be of minor importance in this 
regime. 

We expect that the small discrepancy between the 
measures for /z is caused by the lensing effect of the con- 
densate, which slightly alters the density profile. If the 
imaging lens is aligned incorrectly, more lensing is ob- 
served and the discrepancy between the chemical poten- 
tial based on the phase and sizes increases. 
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D. Observation of interactions between the 
thermal cloud and the BEC 

The effect of interactions between the thermal cloud 
and the BEC becomes important if a significant ther- 
mal fraction is present. As pointed out in Ref. [19], the 
thermodynamic behavior of the cloud is fixed by two pa- 
rameters: the reduced temperature t — T/T^ and the 
ratio rj given by 



TF 



N 



1/6' 



2/5 



(34) 



Here r° = fujj{N / (,{2,)^/'^ denotes the transition temper- 
ature in the absence of interactions and a ~ yjhj (muj) is 
the harmonic oscillator length. 

Under our typical conditions rj — 0.3, which is less 
than its value in other experiments studying the interac- 
tions [IHJEQ]. There it is shown that interactions between 
the atoms shift the transition temperature Tc downward. 
Ref. [H] measured the effects of the interactions by spa- 
tially separating the thermal cloud from the BEC using 
Bragg spectroscopy to avoid the need to incorporate in- 
teractions in the description of the measured density pro- 
files. The separation of both clouds depends sensitively 
on the applied Bragg pulse, and introduces an extra un- 
certainty. Furthermore, the accuracy is limited by the 
absence of an exact theory describing the expansion of 
the BEC. 

Using PCI these interaction can be measured directly 
and more accurately than in these previous studies, al- 
though r] is smaller. The measurements are conducted 
as follows. Using evaporative cooling we obtain a cold 
cloud of atoms at a temperature T below under the 
same experimental conditions as the measurement shown 
in Fig. [5] A phase contrast image of this elongated cloud 
is taken with each cloud cooled to different temperatures. 
Each image is fitted to all three 2D models described in 
§ |IV A| and yield a temperature and chemical potential. 
We find the best fit can be made, determined by the 
smallest sum-of-squarcs, using the Popov model and the 
semi-ideal model. The noninteracting model gives infe- 
rior results and cannot be used to accurately describe 
the measured profiles. This has already been noticed in 
previous work [HI EQ] . 

In this experiment we can also see the effect of the 
mean-field potential of the thermal cloud on the con- 
densed atoms although less prominent than the modifi- 
cation of the thermal distribution. This effect is already 
hinted at in Ref. ^18) . The density distribution of the 
thermal cloud causes an effectively larger potential for 
the condensed atoms. The effect of the resulting distri- 
bution is analogous to increasing the trap frequency: the 
cloud gets smaller and denser. Since we can measure 
both density and size separate and with a high accuracy, 
we measure the effect of the interactions by comparing 
the chemical potential based on the axial size and the 
accumulated phase. The axial size is favored above the 
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FIG. 6: The ratio /iax/^c/) of the chemical potential based on 
the axial size i?ax and phase <j) as function of the temperature. 
Figure (a) corresponds to = 0.3 and (b) corresponds to 
rj — 0.25. In the second series is reduced by reducing the 
number of atoms. 



radial size for two reasons: lensing effects in this direc- 
tion are negligible and the size is large compared to the 
resolution of the imaging setup. 

From the fit results we determine the axial size and 
phase of the condensed part of the cloud and calculate 
its chemical potential based only on the accumulated 
phase, yielding /i^ or the axial size, yielding /^ax- The 
ratio /^ax//^0 is plotted against the temperature T and 
shown in Fig. |6] We find the chemical potential based on 
the axial size for temperatures close to Tc to be approxi- 
mately 15% too small compared to the chemical potential 
based on the accumulated phase. As the temperature de- 
creases, this effect gets smaller until at low temperatures 
(low thermal densities) we find correspondence between 
both measures within a few percent, consistent with the 
results found in the previous section. Since our standard 
fitting procedure uses /x and T to fully describe the den- 
sity profile, we find the sum-of-squares to be 10% smaller 
at temperatures where a significant thermal fraction is 
present (condensate fractions up to 80%) for the fit using 
the Popov model compared to a fit using the semi-ideal 
model. 

Since this shows that the Popov model yields the most 
consistent results all measurements in the remainder of 
this paper are analyzed using the Popov model. 



E. Accuracy and reproducibility 

Repeated measurements are done in the regime where 
three-body losses limit the density of condensed atoms. 
The total number of condensed atoms A^c is therefore 
approximately constant and the results are used to de- 
termine the sensitivity as well as the reproducibility of 
the PCI method. The measurement series takes over 
two hours and during this time the temperature of the 
cloud changes from T = 0.573 ± 0.004/iK to T = 0.387 ± 
0.004/iK due to the apparatus heating up. However, in 
the density limited regime the number condensed atoms 
is roughly constant in this temperature range. The mea- 
sured number of condensed atoms is found to change 
less than 3% during the two measurement time on a 
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FIG. 7: 

Nc and the condensate fraction as a function of the tempera- 
ture T. Each point is the resuh of a fit to the Popov model on 
a single measurement. The statistical uncertainty is smaller 
than the size of the symbol. 



cloud containing N^. = (250 ± 10) • 10® condensed atoms. 
For measurements shot back-to-back the number of con- 
densed atoms changes approximately 1% . Since we as- 
sume the largest contribution is caused by the difference 
in the environmental conditions, this result sets the up- 
per limit of 1% for the shot-to-shot reproducibility of the 
PCI method using the two parameter fit. 

Our fit procedure, in which /i and T are determined 
by both phase and size simultaneously using the Popov 
model, yields statistical errors below one percent, smaller 
than the estimated systematic error. The uncertainty in 
the number of condensed atoms, which scales as iVc ~ 
/i^/^ , is therefore estimated from the discrepancy between 
/.tax and which is roughly 3% and yields an uncer- 
tainty in the number of condensed atoms of 5%. 

In an experiment we set up to determine the accuracy 
of the method on a single shot for various temperatures, 
we produce clouds at various temperatures below Tc- The 
resulting number of thermal and condensed atoms, as 
well as the condensate fraction determined on the fitted 
/i and T is shown in Fig. [7] where each point corresponds 
to a single measurement. The fluctuations in these results 
are only a few percent for both the number of thermal 
and condensed atoms and this measurement shows the 
accuracy of the PCI method over the full temperature 
range. 



F. Transition temperature for condensation 

Generating clouds above as well as below the transition 
temperature allows us to determine the temperature at 
which condensation sets in. The Popov model yields /i 
and T in both regions, where fi becomes negative for 



temperatures above T°. Note that our implementation 
of the Popov model reduces to the Hartree-Fock model 
for T > Tc [16 . 

For T— T° ^ the chemical potential as a function of 
temperature for a noninteracting Bose gas can be written 
as M 



(35) 



where is the Riemann-Zeta function. Equation (35 1 can 
be used to estimate T°. Note that the repulsive interac- 
tions lower the central density and therefore reduce the 
transition temperature. Therefore, condensation does 
not set in at , where ^ — 0, but for /i = 2?7o^gx- 

In the experiment we generate clouds at different tem- 
peratures above as well as below the transition tempera- 
tures and derive T and fj, from a fit to the Popov model. 
The results presented in Fig. |8] show positive values of 
fi for all temperatures, although no BEC is observed in 
the five points with the highest temperatures. The seven 
data points with t he highest temperature are used to 
make a fit to Eq. ( [ssj ) with being the only free pa- 
rameter. This yields r° = 0.998 ± 0.007/iK. The ther- 
mal density is determined by averaging the two measure- 
ments which are the closest to the temperature where 
a condensate is formed. Both rip^ and T? are used to 



find the temperature in Eq. (35), where /i — 2C/o"ox — 
h{1.65 ± 0.05kHz). This yields Tc = 0.972 ± O.OOS/fK. 

The measured shift due to the interactions is {T^ — 
rc)/rc = 0.027 ± 0.001, which agrees well with the the- 
oretical prediction based on the Hartree-Fock mean-field 
theory, ATJT^ w -l.SSa/aN^^^ « 0.025 [B]. 

In order to confirm the systematic errors are well es- 
timated, we compare the value of Tc to the highest tem- 
perature where we observe a bimodal distribution in the 
image, both by fitting the cloud and by examining line 
profiles through the center of the cloud. The highest 
temperature where a BEC is observed is determined at 
T = 0.963±0.005/iK, where the lowest temperature with- 
out a BEC yields T = 0.975 ± 0.003/iK and the deter- 
mined value of Tc indeed lies in between these two tem- 
peratures. The consistency of these numbers suggest the 
systematic errors are well estimated. 



VI. CONCLUSION 

In conclusion, this paper describes the accurate deter- 
mination of the density distribution of condensates at 
finite temperatures using PCI. A detailed description is 
given of the amount of phase a probe beam accumulates 
when a cloud of atoms is passed. The resulting images 
show a periodic variation of the intensity as a function of 
the accumulated phase and the resulting high dynamic 
range is used to measure both the density distribution 
of a thermal cloud and a condensate simultaneously at 
high accuracy. Since all measurements are done in situ, 
there is no need to describe the expansion of the cloud. 
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FIG. 8: The chemical potential a cold cloud derived from the 
fit to the Popov-model as a function of temperature, both 
above and below Tc. The line is the result of a fit to Eq. ( |35[ ) 
on the seven data points where fi/h < 2.5kHz. The data 
point accentuated by the box marks the observation of a small 
BEG. The point accentuated by the cross is the point with 
the lowest temperature, in which no bimodal distribution is 
observed. The dashed line indicates /i = 2Uoncx, the expected 
value for condensation to set in, calculated using the average 
thermal density of the two accentuated points. 

The size of the condensate, also a measure for its density, 
can therefore be determined accurately as well. The size 
of the thermal cloud, the measure for the temperature 
of the cloud, is also determined from the same image. 
Only two relevant parameters, jj, and T, are needed to 
describe the measured clouds given the harmonic confine- 
ment as /i and T determine both the peak density and 



size of both components of the cloud. This procedure 
is used to discriminate between three models describing 
the equilibrium properties of trapped clouds at nonzero 
temperatures, each model incorporating the interactions 
between thermal atoms and condensed atoms differently. 

We resolve the effect of the interaction on the density 
distribution of thermal cloud as well as the smaller ef- 
fect the interaction has on the density distribution of the 
condensate. The model accounting for these effects gives 
consistent results over the complete temperature range. 
The imaging scheme is used to determine the number of 
atoms within one percent from shot-to-shot at our typical 
number of condensed atoms between 2 • 10^ and 3 • 10^. 
Systematic errors are estimated by comparing the dif- 
ferent measures for the chemical potential and we find 
the uncertainty in the number of number of condensed 
atoms to be roughly five percent. PCI can be used both 
below and above the transition temperature and we used 
it to determine and derived the shift of the transition 
temperature due to the interactions. The agreement with 
the theoretical value of this shift indicates the systematic 
errors are smaller than the uncertainty in the determina- 
tion of T. 
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